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Abstract. We present a comparison between the bosonization results for quantum quenches and exact 
diagonalizations in microscopic models of interacting spinless fermions in a one-dimensional lattice. The 
numerical analysis of the long-time averages shows that density-density correlations at small momenta tend 
to a non-zero limit, mimicking a thermal behavior. These results are at variance with the bosonization 
approach, which predicts the presence of long-wavelength critical properties in the long-time evolution. 
By contrast, the numerical results for finite momenta suggest that the singularities at 2/cf in the density- 
density correlations and at fc_F in the momentum distribution are preserved during the time evolution. 
The presence of an interaction term that breaks integrability flattens out all singularities, suggesting that 
the time evolution of one-dimensional lattice models after a quantum quench may differ from that of the 
Luttinger model. 
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A common protocol to study out-of-equilibrium prob- 
lems is called quantum quench and consists in preparing 
the system in the ground state of a given Hamiltonian 
and then suddenly let it evolve under the action of a new 
Hamiltonian. Since the evolution is unitary, the energy 
stored into the initial state is conserved during the dy- 
namics. The interest in these classes of non-equilibrium 
problems relies on both the dynamics itself [5] and the 
long-time properties, including the highly debated issue of 
thermalization [7] . Quantum quenches have been the sub- 
ject of vast literature focusing on different systems [5]. In 
particular, one-dimensional (ID) models have been largely 
explored [5] . From the theoretical point of view, there are 
many different analytical or numerical methods that may 
give useful insights into the dynamical properties and the 
nature of the steady state (if any), inquiring the possibil- 
ity to reach thermalization |9I10I11|I2] . Intcgrablc models 
require a separate analysis |13j . because a complete ther- 
malization is not expected, due to the existence of a ex- 
tensive number of conserved quantities. In this regard, it 
has been suggested that the steady state of an integrable 
system of hard-core bosons may be described by a gener- 
alized Gibbs ensemble (GGE) that maximizes the entropy 
with all possible constraints imposed by the existence of 
(infinite) integrals of motion [14] . 

At the roots of many-body theory in condensed matter 
theory, the Luttinger model was introduced to describe a 
system of interacting fermions in ID; the approximation 
of considering fully linear dispersions made it possible an 
exact solution in terms of bosonic variables, e.g., fcrmionic 



quenches in one-dimensional gapless systems 

densities [15116] (hence the name bosonization) . Later, the 
asymptotic forms of one- and two-particle correlations in 
equilibrium were obtained by Luther and Peschel [17] and 
Haldane [18119] proposed that this model may generically 
describe the low-energy properties of a broad class of sys- 
tems in ID, now known as Tomonaga-Luttingcr (TL) liq- 
uids [20]. 

Recently, the bosonization approach for the Luttinger 
model has been used to compute the time evolution af- 
ter a quantum quench [5]. According to these results, real 
space correlation functions evolve towards a steady state. 
Remarkably, the spatial decay of correlations is always 
governed by critical exponents, which however are dif- 
ferent from the ground state ones. On the other hand, 
generic quantum models in ID are expected to give rise 
to thermalization, where all singularities are washed out. 
It is therefore important to test the bosonization results 
against numerical simulations of lattice models in order to 
better understand the non-equilibrium dynamics of quan- 
tum one-dimensional models. Recently, a numerical in- 
vestigation of a one-dimensional lattice model of spin- 
less fermions has been performed by using density-matrix 
renormalization group technique [21] . The analysis of the 
short-time evolution of few observables suggests that the 
bosonization predictions arc indeed verified also in lattice 
models. By contrast, a renormalization-group approach of 
the bosonized Hamiltonian in presence of a periodic po- 
tential has suggested that temperature effects and dissi- 
pation among bosonic modes may be generated, removing 
all singularities [22] . In this paper, we aim at perform- 
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ing a detailed comparison between the TL approach and assumed that well defined values for the forward scatter- 
exact calculations (by Lanczos diagonalizations) on a mi- ing (small g's) and backscattering (q rj ±2kp) are given, 
croscopic model of interacting fcrmions on a ID lattice. The basic idea is to linearize the dispersion relation 
Although the latter approach is limited to relatively small near the Fermi energy, i.e., ~ uf(|/c| — kp). It is con- 
system sizes, we obtain a clear and unambiguous evidence venient to assume that this linear dispersion extends for 
that a TL analysis misses important aspects of the long- all k G [—00,00]. In this case, the particle belonging to 
time behavior. either branch, denoted by left or right, are distinguishable 
The paper is organized as follow: in sections[5]and[3J we and will be considered as two different species of fermions 
review the basic steps of the bosonization technique and (Luttinger model). [TS] 

the results for the time-dependent correlation functions; in The Hamiltonian (TTJ) can be written in terms of two 

section^ we present the lattice models and the numerical density operators (for left and right particles): 

method based upon Lanczos diagonalization; in section [SI / N t 1 \ t 1 »\ 

PL{q) = }_^ C k C k+q PB\q) = 2^ C k C k+q- (4) 

we make some preliminary considerations on time evolu- 



fc<0 fc>o 
Indeed, up to a constant, we obtain: 



tion on finite sizes and time averages; in section[6j we show 

the results; in sectional we draw our conclusions. Finally, Ua = ~L~ ^ b^W"?) + Ph(q)ph{-q)\ , (5) 



in the appendixEl we show that the exact time evolution ^ = V_ ^ + [()R (_ q) + ph{ _ q)] ( (6) 

of the TL model can be described by a density-density 



V_ 

2L 

q>0 



Jastrow wave function. 



2 Bosonization 



The TL model is obtained from a generic Hamiltonian of 
interacting spinless fcrmions in a ID lattice: 



where V = V(0) — V(2kp). In the standard approach for 
bosonization, new boson operators are considered: 

/ 27T , . , + / 2tt . . , . 

q = \/^L PL ^ q) b ^ V^ Pi(<Z) ' ( ) 

such that a and b obey to the usual canonical bosonic 

commutation relations. In terms of these operators the 

Hamiltonian becomes: 



H — H.Q + Hint, (1) 

^o = J>4^, (2) Wo = «F5>(aJa g + &t& (r ) (9) 

k i>° 

Hint = ^2 V{q)clcl,ck>- q ck +qi (3) "Hint = ^ q K a « + b % + a i b i + b l a l) ■ ( 10 ) 

k.k'.q^O q>0 

. j. . it-,, Therefore, the full Hamiltonian may be easily diagonalized 

where represents the lermiomc dispersion and L is the 

™ .„ „ r,i-, • it i \ by using a q- independent Bogoliubov transformation: 

size ol the lattice. I he specific form ol the interaction V(q) 

depends upon the details of the model. In general, it is a q — ua q + vb^ q , (II) 
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j3 q = va\ + ub„ 



(12) and therefore 



with u 2 — v 2 = 1; therefore, we can take u = cosh(/> and 
v = sinh0, with tanh(20) = V/(2irvp + V). We mention 
that u and v can be expressed in terms of the standard 
Luttingcr parameter K, namely u = (1 + K)/(2VK) and 



-mt 



|0) a6 oc exp 



u ^ 

q>0 



|0)a/ 



cx exp 



L <?>0 



i) 4$ 



\0)ab 



(17) 



In appendix [21 we show that both the exact ground 
v = (1-K)/(2VK), with K = (u-v) 2 = cosh2</>-sinh20. state flgj an d the time evolution of Eq. (JT2J may be 
After the Bogoliubov transformation, the TL Hamil- rewritten in term of density-density Jastrow wave func- 



tonian becomes 



<?>0 



(13) 



where: 



v F V 



q = c s q 



is the excitation energy for the g-mode. 

Finally, let us indicate by |0) a f, and \Q) a p the vacuum 
states of (a, b) and (a, j3) bosons, which coincide with 
the ground states for the non-interacting [Eq. ©] and 
the interacting [Eq. (|13|) ] systems, respectively. Then, the 
ground state of the full Hamiltonian can be written as: 



\0) a p oc exp 



u ^ 

q>a 



tflt 

9^9 



|0> o6 



(15) 



Let us now consider the time evolution of the non- 
interacting state \0) a b by using the interacting Hamilto- 
nian "H, corresponding to a quantum quench from V = 
to a finite value of the interaction strength. Interestingly, 
the time evolution has a very simple and instructive form, 
which will be used in the following Section. Indeed, from 
Eq. (TPo!) we can write 



\0)ab oc exp 



q>0 



|0>c 



(16) 



tions, [23] which are commonly used to describe ground- 
state properties of correlated systems, [24125] and more 
recently also out of equilibrium dynamics. [26] 



( 14 ) 3 Correlation functions 



Here, we give a brief overview of the time-dependent cor- 
relation functions; our results agree with those of Ref. [5] 
for the Luttinger model. Let us suppose that at t = the 
system is set in the non- interacting ground state |0) Q t>, 
while for t > it evolves according to the interacting 
Hamiltonian %: 



(18) 



Then, the bosonization technique allows one to directly 
compute the density-density correlation function: 



N q (t) = i<*(t)|n,n_,I#(*)), 



(19) 



where n± q — p(±q) = ph(±q) + Pi?(±<?)- Indeed, by using 
Eqs. 0, flgj), flU), ([12} and ([17]). after some straightfor- 
ward algebra, we obtain the simple expression: 



Nq(t) = ^— [A + B cos(2e g *)] 

Z7T 



(20) 
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where 4 Models and numerical methods 



A = (u 2 + v 2 )(u~v) 2 = ^-^, (21) 



B = l-A. (22) 

This result shows that N q (t) docs not relax to a finite 
limit but instead oscillates with a single frequency that 



Let us now define the microscopic model that will be used 
to make comparisons with the bosonization predictions: 
a system of spinlcss fcrmions interacting via a repulsive 
short-range potential: 



equals twice the value of the excitation energy e q = c s q. By % = _ J ^ c | c . + Hc + y £ n . n . +V > £ n . n . ( 26 ) 



performing a Fourier transform (that takes into account 



an ultraviolet cutoff), it is possible to obtain the q ~ where and indicate nearest-neighbor and next- 

(spatially monotonic) contribution of the density-density nearest- neighbor sites, respectively, cj (a) creates (de- 
correlations: stroys) a fermion on the site i, and rij = c\ci is the fermion 

density. The number of sites and fermions are denoted 

. (23) 



*r(t)=n a -JL 



2A B B 



(r + 2c s t) 2 (r - 2c s t) 
which shows that, for each fixed distance r, lim^oo N r (t) 



by L and N, respectively, so that the fermion density is 
n = N / L. Periodic or antiperiodic boundary conditions 
are considered for odd or even N, respectively. In the 

converges to a finite value. Notice that the asymptotic 

following, we take J = 1, as energy scale. Most of the 

behavior of N r (t) coincides with the Fourier transform of 

calculations are done in presence of nearest-neighbor in- 

the time average of N q (t), which equals Aq/2-K and differs 

teraction V only (which corresponds to an integrable sys- 

from the known ground state behavior Kq/2iT. 

tern), similar results are also obtained including a finite 

For completeness, wc also report the time evolution of 

next-nearest-neighbor V (that breaks the integrability). 

the q ~ 2kp singular contribution to the density correla- 

Notice that the system we investigate is equivalent to a 

tions: 

model of hard-core bosons with Hamiltonian (|26[) and pe- 

/ a 2j-2 \ 2uv(u — v) 2 v 

N (t) oc ( s | cos(ZKp r) riodic boundary conditions. 

r ^ ' ^| r 2_4 c 2^2|y r 2(u 2 +v 2 )(u-v) 2 ' <■ ' 

The quantum quench consists in taking an initial wave 

and the q ~ k p singularity of the one-body density matrix: . . 

function Itf'(O)), which is the ground state of Eq. (|26[) with 

, , .. + . , ,. f\r 2 -4c 2 t 2 \\ 2u2v2 cos(k F r) . , V = V t (and V = V{), and letting it evolve under the 

m\4c mt)) oc ^ / ' j ^^y. (25) 

s same Hamiltonian with V = Vf (and V = Vj). By using 

At fixed r, for t — > oo the first factor in these expressions the Lanczos method, it is possible to perform the exact 

tends to unity, showing that the correlation functions in time evolution of any initial state. Indeed, the full time 

real space tend again to a finite limit for large times. interval can be split in small steps At and the time evo- 
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lution can be evaluated recursively: where exact results can be obtained in the thermody- 

namic limit. In particular, we consider the case where the 

\V(t + At)) = e- mAt \V(t)}. (27) 

density of spinlcss fcrmions is half filled (i.e., n = 1/2) 

Each small-time evolution can be computed by a trim- an d the initial state has a fermion every two sites, i.e., 

cated Taylor expansion: |#(q)) = | . . . , 1, 0, 1, 0, 1, 0, . . .). Then, the time evolution 

(-iAt) k , is done by the Hamiltonian (EBJ with V = V = 0. In this 

\#(t + M))~£L-^H k Mt)), (28) ^ 



fe=0 

where the cut-off k c is chosen as to preserve energy con 
servation to the desired numerical accuracy. 



case, density-density correlation functions 



N r (t) = (&(t)\n r n \&(t)) - (V(t)\n r \#(t))(*(t)\no\V{t)) 

(30) 

can be computed analytically directly in the thermody- 
namic limit. This correlation function is translationally 
From what we have described in the previous section, by invariant and its Fourier transform is given by: 



5 Preliminary considerations 



using the Lanczos technique, it is in principle possible to i 

N q (t) = -[l-M8tsm(q/2))}, (31) 

compute the time-evolution of any observable: 

where Jo (a;) is the Bessel function of order zero. In Fig. [I] 

0(t) = (*(t)\Omt)}. (29) 

we compare Eq. ([3"TT) with the case of L = 32 for q = 
However, we must emphasize that, on finite systems, re- tt/16. Although the results for L = 32 clearly deviate from 
currence effects are always present and, at long times, the analytical ones for t > 35, their time average is fully 
the dynamics suffers from size effects. Indeed, the ini- meaningful and gives an excellent approximation of the 
tial state inevitably contains excitations that, under the exact outcome. Therefore, in the following we will consider 
time evolution, may "travel" all along the chain. In mod- time averages over relatively large times (up to t = 100) 
els where the full spectrum is described by non-interacting to obtain an estimation of the long-time behavior, 
quasi-particles, interference effects appear after the time 
T ~ L/v, where v is the velocity of the elementary exci- 

6 Numerical results 

tations. However, for generic models, this recurrence time 

may be much larger and relies on the energy differences We begin by considering the case of V = 0. In this case, 
of the (finite-size) many-body spectrum. Nevertheless, al- the low-energy properties of Eq. (|26[) are well described by 
though for large times 0(t) may suffer from size effects, a (gapless) TL liquid, except for n = 1/2 and V > 2, where 
its average over long times is fully meaningful. In order a (gapped) charge-density-wave insulator is obtained. |20j 
to give some support on this claim, we consider a case Since we want to compare the numerical calculations with 



E. Coira, F. Becca, and A. Parola: Quantum 




t 



Fig. 1. (Color online) Time evolution of the density- 
density correlation N q (i) of the inhomogeneous state 
| . . . , 1, 0, 1, 0, 1, 0, . . .) that is evolved with the non- 
interacting Hamiltonian of Eq. ([26]) . Results with L = 32 
(solid black curve) are compared with the thermodynamic 
limit (dotted black curve). The time averages are also re- 
ported (blue dashed lines). 

the TL theory, in the following, we take n — 1/4 (similar 
results have been also obtained with other densities, e.g., 
n = 1/3). We will show the case of Vi = (and V{ = 0), 
so that the initial wave function corresponds to the case 
of free spinless fermions. However, we checked that quali- 
tatively similar results are obtained also for finite Vi > 0. 

Some initial information about the time evolution can 
be obtained from the analysis of the overlaps between the 
initial wave function 1^(0)) and the eigenstates of the 
evolving Hamiltonian \<P n ) (with H\^ n ) — E n \<P n )) i.e., 
c„ = (^nl^O)). Indeed, these coefficients play an impor- 
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tant role in the time evolution of a generic observable O: 

(#(t)|0|ff(t)) = e-« E »- E »''> t <*,c n ($ n ,\0\$ n ). (32) 

This expression shows that the temporal evolution of O 
contains all the frequencies u> = E n — E n i corresponding to 
the gaps in the excitation spectrum of the interacting sys- 
tem. Bosonization predicts a linear excitation spectrum at 
low energies that, in the special case of the observable N q 
written in terms of density operators, gives rise to a single 
frequency oscillatory behavior with lu = 2e q , see Eq. (|20p. 
The excitation spectrum of a lattice model is consider- 
ably more complex and curvature effects are expected to 
introduce further frequencies in the power spectrum of 
the time evolution of O. Pure oscillatory behavior would 
correspond to a c„ distribution peaked in a small energy 
interval E n — E n i and to matrix elements (<P n '\0\ ( I> n ) able 
to connect a single pair of excitations. 

In Fig. [21 we show the results for P(E) in two differ- 
ent cases with Vf = 3 and 10, which will be considered 
in the following discussion. The c n distribution is indeed 
confined in a limited energy interval although a significant 
additional peak at high energy appears for strong interac- 
tions. 

Let us now move to the main part of the paper and 
consider the density-density correlations: 

N <® = lEe^'-^TOk^lW (33) 

To compare the exact diagonalizations with bosonization 
results of Eq. (f2T)]) , we have to compute the parameters A 
and B, together with the renormalizcd fermionic disper- 
sion e q = c s q. This can be easily done by calculating the 
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Fig. 2. (Color online) Distribution P(E) of the overlaps Fig 3 (Color onlinc ) Timc evolution of the density- 
squared) between the initial wave function and the eigen- density corrclation Nq ^ at quartcr ming for a quench 
states of the evolving Hamiltonian as a function of the from y. = to Vf = 3 and q = 2tt/32 (the system size is 
energy of the eigenstates E. Energies have been shifted, L = 32 ). The timc average (dashed line) and the bosoniza- 
namely AE = E - E i (where E t is the energy stored into tion results ( dotted line ) are also snown . 
the initial state), so that AE = corresponds to the av- 
erage value of the distribution. Lattice sizes are L = 20 
(dashed line) and L = 32 (solid line). A broadening of the 
finite size 5-like peaks has been introduced. The value of 
the interaction strength is V — 3 (left panel) and V = 10 
(right panel). 



0.07 
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0.04 

Luttinger parameter K and the velocity c s , which can be 
obtained cither by Bethe Ansatz (for V' = 0) or numeri- 
cally (for V j£ 0). [ID] 

In Figs. El 01E1 andl we report the results for V f = 3 Fig _ 4 _ (Color online) Samc ag ^ Fig Efor q = 4 ^ /32 
and 10 and the two smallest non-zero momenta on the 32- 

site lattice, namely q = 2ir/32 and 47r/32, together with with the bosonization approach. Nevertheless, long-time 

the bosonization prediction of Eq. ([2"D"]) . The time aver- averages in the TL model are always different from those 

age of the oscillating signal stabilizes at a well defined obtained in the lattice model, where the signal shows sig- 

value for sufficiently large times, a feature also shared nificant deviations from periodicity, as a consequence of 
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Fig. 5. (Color online) Same as in Fig. [3] for V/ = 10. 
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Fig. 6. (Color online) Same as in FigJ5]for q = 4ir/32. 

the coupling among excitation modes. However, for small 
quenches (e.g., Vt = 3) and small momenta (e.g., q — 
27r/32), after an initial transient, a stable oscillatory be- 
havior dominated by a single frequency is observed. In 
this case, the observed frequency agrees very well with 
the TL predictions. For larger values of the momenta and 
for larger quenches, the signal acquires a more complex 
periodicity 



We remark that both the amplitude and the average 
value of the numerical results considerably differ from the 
bosonization predictions. The discrepancy increases for 
larger quenches, where, even for the smallest momenta 
available in the Lanczos diagonalizations, the dynamical 
signal contains more than one frequency and the TL pre- 
dictions become less and less accurate. For example, for 
Vf = 10 and q = 2ir/32, the lattice model shows at least 
two relevant frequencies (see Fig.[S]), the largest one being 
very close to the TL result. Also the discrepancy between 
the average values grows when increasing the final value 
of the interaction strength. 

Although the signal may have a very strong depen- 
dence on the system size for large times, the average value 
is quite stable, showing that the long-time properties may 
be safely extracted from our finite-size calculations, see 
Fig. [7] Moreover, there is some evidence that fluctuations 
around the average value decrease by increasing the cluster 
size, suggesting the possibility that in the thermodynamic 
limit (and long times) the signal experiences a complete 
damping towards its average value. 

From these results, it is clear that bosonization misses 
relevant aspects of the dynamics after a quantum quench. 
More surprisingly, the long-time average of the density- 
density correlations are qualitatively different from the bosoniza- 
tion predictions for q — > 0. Indeed, in the TL theory, a lin- 
ear behavior N q = A/(2n)q (where the overbar indicates 
the long-time average) is found from Eq. ([2"U]) . while, the 
numerical results clearly indicate a finite limit N q = const 
as q — > 0, see Figs. [5] and GO The slope of N q predicted by 
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Fig. 7. (Color online) Time evolution of the density- _. _ ,„ . N _. . , ., . 

v Fig. 8. (Color online) lime average of density-density 

density correlation N a (t) at quarter filling for a quench . Ar/ , Nr lr Tr «, T ^ or 

y correlation -/V g (t) for a quench from Vi = to Vf — 3 for 

from V t = to V f = 3 and g = 2tt/32 for L = 16 (dashed £ = 2Q (squares) ^ £ = 32 (drdes) dashed ^ 

lines) and L = 32 (solid lines). The two averages, eval- 
uated after a long time interval (i.e., t = 100), are also 
reported (horizontal lines). 



shows the behavior expected from bosonization at small 
momenta. 



bosonization is also shown in the figures and remarkably 
describes the behavior of correlations in an intermediate 
range of wave-vectors. It appears that although the TL 
model does not capture an important qualitative feature 
of the long wave-length behavior of correlations, it is still 
able to correctly reproduce the physics of the asymptotic 
state when the size of the system is not exceedingly large. 
We will see that similar features also appear in other ob- 
servables, like the momentum distribution. 

Although Lanczos diagonalizations are performed on 
finite lattices, these conclusions are robust against finite- 
size effects, as proved by the nice collapse of the numerical 
data on L = 20 and 32 clusters, see Figs. [8] and [9] On the 
one hand, a non-vanishing long-wavelength limit of N q is 
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Fig. 9. (Color online) Same as in Fig. [8] for Vt = 10. 
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reminiscent of a finite temperature behavior, where sin- 
gularities are washed out by thermal fluctuations. [37] On 
the other hand, our finite-size calculations provide some 
evidence for the existence of a cusp at q = 2kp. Unfor- 
tunately, due to the intrinsic statistical error induced by 
the time average procedure, the (critical) exponent related 
to the q = 2k p singularity cannot be accurately deter- 
mined from Lanczos data. The possibility that the cusps 
are eventually rounded in the thermodynamic limit cannot 
be ruled out either. In any case, we find that N q cannot 
be suitably fitted by using a single effective temperature 
(as expected, since the model is intcgrable). 

The failure of the bosonization approach may be due 
to the presence of a perfectly linear fermionic dispersion 
(up to infinity or to a given cutoff) in the TL model. While 
this approximation is known to give the correct low-energy 
behavior for static properties, |28I29| it breaks down when 
considering the real time evolution for long times, where 
the initial state may contain high-energy excitations, not 
adequately represented within the TL model. Moreover, 
bosonization describes uncoupled modes, which do not in- 
teract and, therefore, does not include the effects due to 
dephasing. In this regard, band curvature and finite band- 
width effects are expected to play a crucial role beyond the 
simple TL results. 

Let us now consider the momentum distribution: 



n k . 



(*) ^E^^'^l^l^))' ( 34 ) 




The numerical results for the time average are shown 
in Fig. [TU] for Vf = 3 and 10. On any finite-size system, 
the momentum distribution shows a jump at hp, which 



Fig. 10. (Color online) Time average of the momentum 
distribution for L = 20 (squares) and L = 32 (circles) for 
Vf = 3 (bottom panel) and Vf = 10 (upper panel). 

is clearly visible in our numerical results. In the thermo- 
dynamic limit, bosonization predicts the occurrence of a 
non-analytic behavior in n^: An^ oc \k — kp\ a , with a > 0. 
The exponent a could be extracted from a finite-size scal- 
ing analysis, in order to determine the presence (for a < 1) 
or the absence (for a > 1) of a singularity. 

We note that the available sizes are sufficient to ob- 
tain a very accurate determination of the exponent for the 
ground state, where the exact values can be computed by 
Bethc Ansatz. On the contrary, a similar fit for the time 
averaged results appears to be more problematic, as shown 
in Fig. 1111 Indeed, a pure power-law fit of the data appears 
to be less accurate and the resulting exponent a increases 
when the fit is limited to the largest sizes. Nevertheless, 
Fig. [TT] shows that the bosonization approach provides a 
value for a roughly consistent with the numerical analysis, 
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Fig. 11. (Color online) Size scaling of the jump at q = kp 
in the time average of the momentum distribution. Results 
for Vf = 3 (bottom panel) and Vf = 10 (upper panel) are 
reported (full circles). The results for the ground state 
(squares) are also reported for comparison. The slope of 
the lines corresponds to the critical exponents predicted 
by bosonization. 

in agreement with very recent density-matrix rcnormaliza- 
tion group calculations. |21j However, we remark that a 
complete smoothing of the curve at larger sizes cannot be 
excluded by the Lanczos data. 

Finally, we consider the effects of a finite Vf in the 
final Hamiltonian. When sufficiently large, this term is 
expected to break the integrability conditions. In this case, 
thermalization should be expected on general grounds and 
we are in the position to verify whether some evidence 
for that is already visible on the lattice sizes studied by 
Lanczos diagonalization. In Figs. [T^] and 1131 we report 
the results for N q and rife, respectively. Three different 
values of the final interaction strengths are reported. In 
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Fig. 12. (Color online) Time average of density-density 
correlation N q (i) for a quench from Vi = to Vf = 10 and 
V* = 3 (a), V> = 15 and V/ = 4.5 (b), and V f = 20 and 
Vf —6 (c). Data for L = 20 (squares) and i = 32 (circles) 
are shown. Thermal values corresponding to an effective 
temperature that gives the correct internal energy are also 
shown (empty circles). 

all cases, N q = const for q — >• 0, while a sizable peak is 
still present at q = 2fcp. Therefore, the qualitative picture 
of the previous integrable model holds also in this case. 

We also notice that, for small quenches (e.g., Vf = 10 
and Vf = 3), both the density-density correlations and 
the momentum distribution may be nicely fitted by as- 
suming a single effective temperature within the canon- 
ical ensemble, see Figs. [T2] and [T21 Indeed, by choosing 
the temperature that reproduces the value of the internal 
energy, we are able to obtain a quite satisfactory repre- 
sentation for all momenta, at variance with the integrable 
case, Vf = 0. Of course, in order to show that a real ther- 
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Fig. 13. (Color online) Same as in Fig. [12] for the mo- 
mentum distribution rikit). 



malization takes place, one must show that all correlation 
functions are described by a single effective temperature. 
In this respect, the best way is to consider the properties 
of the (reduced) density matrix, as suggested by Poilblanc 
in a recent work. [30] This task is beyond the scope of the 
present paper, which is centered around the comparison 
between the TL model and a microscopic model on the 
lattice. 



In conclusion, we reported a direct comparison of the non- 
equilibrium dynamics between the bosonization approach 
and exact diagonalizations for an interacting model of 
spinlcss fcrmions. On the one hand, we showed that the 
bosonization technique does not capture few important as- 
pects of the long-time behavior. In particular, a thermal- 
like behavior of the density-density correlations at small 
momenta is observed in the numerical calculations (i.e., 
N q ~ const). On the other hand, our numerical calcula- 
tions point towards the occurrence of a singularity at the 
Fermi wavevector in the momentum distribution, as pre- 
dicted by bosonization. Our conclusions do not crucially 
depend upon the presence of a next-nearest-neighbor in- 
teraction V that breaks integrability, and show that the 
critical behavior predicted by the bosonization technique 
should be considered with care and may strongly depend 
upon the linearization of the fermionic band. The numer- 
ical analysis for non-integrable models is consistent with 
thermalization in the thermodynamic limit, however, sig- 
nificant size effects hamper the possibility to quantita- 
tively investigate the thermalization dynamics. 

We thank M. Fabrizio, G. Santoro and A. Silva for 
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For larger quenches, an accurate description of these jj MaNEP) 



correlation functions in terms of a single effective temper- 
ature is less accurate (see Figs. [T2l and [T3l and quite sub- 
stantial deviations from the thermal values are observed 

for large momenta in the density-density correlations and In this Appendix, we show that the ground-state wave 



A The Jastrow wave function 



for small momenta in the momentum distribution. 



function (|15|) of the Tomonaga-Luttingcr Hamiltonian may 
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be rewritten as a Jastrow term acting on the non-interacting 
ground state \0) a b, i.e., in the form: [3T] 



\0) a/3 oc exp 



q>0 



jqi tqf L—q 



\0)ab 



(35) 



where the density operators are defined by n± q — p(±q) = 
p^(±q) + pji{±q). Most importantly, we also show that 
the exact time evolved state (|17[) can be also written as a 
(time-dependent) Jastrow term acting on \0) a b- 

Let us start with the ground state. In order to find 
out the expression of the pseudo-potential w q , on the one 
hand, we write \0) a p by using Eqs. (fTTj) and (|T2l) : 



|0) Q(3 oc exp 



^2( ua l + vb q )(va q + ub\) 

U q>0 



\0)ab- (36) 



On the other hand, we use the fact that: 

n q n- q = ^-{a q a\ + b\b q + a q b q + b\a\). (37) 

Since, in a Tomonaga-Luttinger liquid, different modes la- 
beled by q > are not coupled, we can focus our attention 
on a generic mode, dropping the label q. We want to see 
whether, for a suitable choice of the amplitude w q , the 
expression in Eq. (|3"5"j) coincides with that of Eq. (131)1) . A 
solution exists if it is possible to define an amplitude / 
such that the following equality holds: 

e /(a t a+k 1 'h+ah+h t a t )|Q^ ^ £ x(aa f +b f b+Dab+±tf(J) |q\ 

(38) 

for all choices of x and v. In the case of interest, due to 
Eq. ([31)f we have x = —v 2 , v — v/u, while Eqs. (l35ll . 
and (|5T)) give / = w q q/(2n). Via some lengthy algebra it 
is possible to show that Eq. (|3"8"|) is satisfied by the choice: 



q(u — v) ' 



(39) 



in agreement with the results of Ref. [3T]. The explicit 
form of the Jastrow pseudo-potential shows a long-range 
(density-density) repulsion, with a logarithmic decay in 
real space: 



|0) Q(3 oc exp 



2tt 



/,/ \ n q n —q 

q(u — v) 

q>0 HV > 



\0)ab- (40) 



Let us now consider the time evolution and show that 
the time evolved wave function can be written as a time- 
dependent Jastrow term applied to the non-interacting 
state. This can be easily proved by considering Eq. (|38|l 
with x = v 2 (e~ 2teqt — 1), v = v/u, and / = w q (i)q/(2ir). 
The final result reads as: 



-mt 



\0)ab oc exp 



j^w q (t)i 



q q 



q>0 



|0)«6 



with 



,(*) 



2nv 



q(u — v) 



u(e 



1) 



-2ie„t 



(41) 



(42) 



Therefore, we obtain the important result that the ex- 
act time evolution of |0) a 6 under the action of H can be 
written as a Jastrow wave function, with a complex time- 
dependent pseudo-potential w q (t). 
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